Manipulating and probing Majorana fermions using superconducting circuits 
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Majorana fermions are long-sought exotic particles that are their own antiparticles and satisfy 
non-Abelian statistics. Here we utilize superconducting circuits to construct two superconducting- 
qubit arrays where Majorana modes can occur. Moreover, we propose to use four superconducting 
qubits as the smallest system to demonstrate the braiding of Majorana modes and show how the 
states before and after braiding Majoranas can be discriminated. This provides an experimentally 
realizable, relatively simple system for manipulating and probing Majorana fermions. 



Introduction. — Majorana fermions are particles that 
are their own antiparticles and satisfy non-Abelian statis- 
tics. These long-sought particles have recently received 
considerable interest (see, e.g., [l|-[s|)- It has been rec- 
ognized [9] that a relatively easy to engineer system — 
one-dimensional (ID) semiconducting wires on an 5- 
wave superconductor — can realize a nontrivial topolog- 
ical state supporting Majorana fermions. This state is 
characteristic of ID topological superconductors [4], in 
which Majorana modes can occur without requiring the 
presence of vortices in the system. Moreover, it was pro- 
posed [10] to use tunable ID semiconducting wire net- 
works on an 5-wave superconductor to demonstrate the 
non-Abelian statistics of Majorana fermions, because the 
Majoranas in the semiconducting wires can also behave 
like vortices in a p + superconductor jH, In ad- 
dition, it was also recognized [11] that when a Jordan- 
Wigner transformation is performed, a ID quantum Ising 
model is equivalent to a ID topological superconductor, 
and Majorana modes can also occur therein. Neverthe- 
less, less attention has been paid to this quantum Ising 
model than to ID topological superconductors because it 
was often regarded as a toy model. 

In this Letter, we unambiguously show that such 
a toy model can be realized experimentally by using 
superconducting-qubit arrays. Importantly, supercon- 
ducting qubits can behave as controllable artificial atoms 
and tunable interqubit couplings are also achievable (see, 
e.g., pjj). These distinct advantages of superconduct- 
ing qubits make it possible to construct a tunable ID 
quantum Ising model on wire networks, similar to the 
semiconducting wire networks in [10], to demonstrate the 
non-Abelian statistics of Majorana fermions. We pro- 
pose to use four superconducting qubits as the smallest 
circuit to demonstrate the braiding of Majorana modes, 
and show how the states before and after braiding Ma- 
joranas can be discriminated. This should provide an 
experimentally realizable, relatively simple setup to ma- 
nipulate and probe Majorana fermions. 



Majorana fermions in superconducting circuits. — We 
construct two types of superconducting-qubit arrays (see 
Fig. [T]), which can exhibit Majorana modes. 

(i) Charge- qubit array. For the array of charge qubits 
shown in Fig. [TJa), every pair of nearest-neighbor qubits 
are coupled by a large Josephson junction acting as an ef- 
fective inductance. The non- nearest-neighbor qubits can 
also be coupled via these large Josephson junctions, but 
the interactions are negligibly small. Here we assume 
that all charge qubits are identical and that all large junc- 
tions are equal to each other. When leading terms are 
considered, the Hamiltonian of the charge-qubit array 
can be written as 
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with/i = \E^Y,{l-CgVg/e), Ej{^q) = Ejo cos(7r$^/<l>o), 
and the interqubit coupling is given by 
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Here, e'^/Cj) ^ Ej^ in the considered charging 

regime and Lj = <l>o/27r/c, with Ic = 27rEjc/^o and 
being the flux quantum. The eigenstates of the Pauli 
operator are the charge states |0n) and |ln), corre- 
sponding to zero and one extra Cooper pair in the super- 
conducting island of the nth qubit. The Hamiltonian ^ 
provides an analog to the ID quantum Ising model. 

We now consider the case with the fluxes in all charge- 
qubit loops being tuned to = ^^Oi so that Ej{^q) = 
0, and the interqubit couplings reach the maximum 
t = Lj{7tEjo/^o)'^. Using the Jordan- Wigner transfor- 
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for all qubits and also for all couplers. Moreover, the 
plasma frequency Vp of the coupler is much higher than 
the related qubit energy, so as to have the coupler in the 
ground state When the leading terms are included, 
the Hamiltonian of the flux-qubit array can be written 
as 
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FIG. 1: (Color online) Two arrays of superconducting qubits. 
(a) Charge-qubit array: Nearest-neighbor charge qubits Qn 
and Qn+i are coupled by a large Josephson junction with cou- 
pling energy Ejc (shown as a crossed rectangle), (b) Flux- 
qubit array: Nearest-neighbor flux qubits are coupled by a 
coupler consisting of a flux-biased loop that is interrupted 
by two large Josephson junctions (each with coupling energy 
Ejc) and a small Josephson junction with coupling energy 
l3Ejc, where < /3 <^ 1. In (a) and (b), $g is the flux applied 
to each qubit loop, (c) Main components of a charge qubit, 
where a superconducting island (denoted as a solid circle) is 
connected to two Josephson junctions (each with coupling en- 
ergy Ejo <^ Ejc and capacitance Cj) and biased by a voltage 
Vg through a gate capacitance Cg <C Cj. (d) Main compo- 
nents of a flux qubit, where two Josephson junctions, each 
with coupling energy Ej <^ Ejc, connects a symmetric dc 
SQUID biased by a flux 



where the Dirac fermions obey the anticommutation re- 
lation {ctniCLn'} = ^nn' - Introducing Majorana fermions: 



one can rewrite the Hamiltonian ([3]) as 
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where = {7^,7^ } = 26xx'Snn^. 

(ii) Flux-qubit array. Figure [TJb) shows an array of 
flux qubits. Here the small junction in the ordinary flux 
qubit is replaced by a symmetric dc SQUID to increase 
the tunability of the qubit. Also, a coupler consisting of 
three Josephson junctions is used to produce a control- 
lable interqubit coupling between nearest-neighbor flux 
qubits. We assume that the parameters are the same 
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Here £ = /p<l>o(| — /), with Ip being the persistent cur- 
rent of the flux qubit and / = ^g/^o + /s/2, where 
fs = ^s/^O' The eigenstates of the Pauli operator 
are the clockwise and anti-clockwise persistent-current 
states of the nth qubit. The symmetric SQUID pro- 
vides an effective Josephson junction with coupling en- 
ergy aEj^ where a = |cos(7r/s)|. The exact expression 
of /i in Eq. ([6]) cannot be given, but it depends on a; nu- 
merical results [15] and approximate analytical calcula- 
tions [16] showed that /i = when a = 1. The interqubit 
coupling strength reads Q 
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where = 2/3 sin(27r/c)/[l + 2/3 cos(27r/c)], and = 
^c/^Q is the reduced flux applied to the coupler. 

We study the case with / = ^ for all flux qubits, 
so as to have e = 0. The Hamiltonian of the system 
also becomes Eq. (|3j) when applying the Jordan- Wigner 

transformation: an = cr+ n^"l\ ^r^. «i = (^n iXn^i c^m. 
where = ^{cr^ ±zcr^). Finally, the Hamiltonian is de- 
scribed by Eq. (jSj by introducing Majorana fermions in 
Eq. (|4|). Therefore, the resulting Hamiltonians in terms 
of Majorana fermions are the same for both charge- and 
flux-qubit arrays. 

For ^ 00, we can obtain the energy bands of the pe- 
riodic chain by performing a Fourier transform on Hamil- 
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tonian © : 7^ = ^ ^ E/c e'^^'lk ^ where 7_^ 
X = B. The resulting Hamiltonian in the reciprocal 
space reads 
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with D{k) = te^^ + /i. The energy spectrum shows the 
particle- hole symmetric dispersion E{k) = ±\D{k)\ = 
±\te'^^ + m|, which consists of two bands. 

When |r| < 1, where r = /i/t, there are two degenerate 
edge modes with zero energy (i.e., in the middle of the en- 
ergy gap) for a finite chain. These two edge modes can be 

represented by Q = Ci7j^+C27f H hC2Ar-i7^+C2Ar7^, 

with the coefficients determined by fiC2n-i +^C2n+i = 0, 
and tc2n-2 + MC2n = 0, where n = 1,2,..., A', and 
the initial condition is cq = for the left-end edge 
state and C2Ar+i = for the right-end edge state. It 
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can be derived that the left-end edge mode is given by 
Ql = C[^f - ni + r'l^ -■■■ + (-r)^- 
right-end edge mode becomes Qr = C[{—r)^' "7r 
^^7^-2 ~ ^7^-1 + 7^]' where the normahzation factor 
is C = {^^=0 '^'^^)~^^'^ ' particular, when /i = 0, the 
Hamiltonian is reduced to = — Yln=i '^i'^^n^n ~ 1)^ 
where dn = ^(7n+i + ^In) ^ Dirac fermion composed 
of two Major anas at adjoining sites. The edge modes be- 
come two isolated Major ana fermions at the left and right 
ends of the chain: = 7^, and Qr = 7^, which do not 
appear in the Hamiltonian because these modes have zero 
energy |3, [l^. Now define \F) to be the state in which 
all eigenstates of the system with E < are occupied 
and those with E >0 are empty. When the edge modes 
are occupied, = and = J§\F) are two 

degenerate ground states of the system. These two Ma- 
jorana modes can be used to represent the states of a 
topologically-protected logical qubit [10]: |0) = dend\F) 
and |1) = dl^^\0), where dend = ^{li ^ '^In) is a non- 
local Dirac fermion, and (iend|0) = 0. 

Manipulating and probing Majorana fermions. — The 
superconducting-qubit arrays proposed above can be 
used to realize a tunable ID quantum Ising model on 
wire networks, similar to the semiconducting wire net- 
works in [10], to demonstrate the non-Abelian statistics 
of Majorana fermions. In particular, braiding Majoranas 
can be implemented via a T-junction formed by two per- 
pendicular wires [10]. Here we use four superconducting 
qubits, as the smallest size of the system, to form such 
a T-junction [see Fig. [2] (a)], where Ej = {e = 0) for 
all charge (flux) qubits. When the Jordan-Wigner tran- 
formation is performed, this T-junction of four qubits is 
described by 

H = t{ai - a\){al + + t{ay - a\,){al + 
-\-t{a2 — <^2)(^3 + ^3) ~ /i(2ajai — 1) 
—fi{2a\,ay — 1) — fi{2a2a2 — 1) — lJi{2a\a^ — 1) 

= ^(^7f 72'' + t^vl2 + ^72''73'') 

-^(^71^7? + ^vlv + M72''72'' + M73''73''), (9) 
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where qubits are numbered by starting from sites 1 and 
V and ending at site 3. 

For Hamiltonian (|9]), when t = for all pairs of 
nearest-neighbor qubits, no edge states occur and the 
whole T-junction region is non-topological. Starting from 
this phase, we adiabatically vary the parameters of su- 
perconducting qubits to have the horizontal array be- 
come a topological region, i.e., the Hamiltonian of the 
system becomes H = i{t'^i^2 + ^lils) ~ ^/^Ivlv — 
—t{2d\di — l)—t{2d2d2 — l)—/J^{2a\,ai' — l). This creates a 
pair of Majoranas at the two ends of the horizontal array 
[see Fig. El^b)]. Generally, the system takes the super- 
position state of these two degenerate Majorana modes: 
1^) = (a7j^ + 67f where \F) is the state in which ah 
eigenstates of the system with E < are occupied. How- 
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FIG. 2: (Color online) Braiding two unpaired Majorana 
fermions. (a) T-junction formed by qubits 1, 2, 3 and 1^ 
where each qubit is denoted by a rectangular box. Two 
Majoranas related to the same qubit (e.g., ji and ji) can 
be paired by the parameter /x while two Majoranas related 
to adjoining qubits (e.g., and 7^) can be paired by t. 
(b) Topological region for the whole horizontal array, where 
a unpaired Majorana (denoted as a solid circle) occurs at 
each end. (c) Adiabatically tuning /i to nonzero for qubit 

1 and turning off t between qubits 1 and 2 drive the left- 
end Majorana mode (shown in yellow) to the middle qubit. 
(d) Adiabatically tuning /i to zero for qubit l' and turning on 
t between qubits 2 and l' drive the original left-end Majorana 
to the bottom of the T-junction, (e) The right-end Majorana 
(in black) is driven to the middle qubit by adiabatically tun- 
ing /i to nonzero for qubit 3 and turning off t between qubits 

2 and 3. (f) The original right-end Majorana is finally driven 
to the left end by adiabatically tuning /x to a sufficiently large 
value for qubit 2, tuning /i to zero for qubit 1 and turning 
on t between qubits 1 and 2. (g) The Majorana at the bot- 
tom is driven to the middle qubit by adiabatically tuning /a 
to zero for qubit 2 and tuning /i to a sufficiently large value 
for qubit l^ (h) Adiabatically turning off t between qubits 2 
and l\ tuning /x to zero for qubit 3, and turning on t between 
qubits 2 and 3 finally drive the original left-end Majorana to 
the right end. This accomplishes the ant i- clockwise braiding 
of two Majorana fermions. 



ever, while reaching the state in Fig. [2fb), if /i for qubits 
1, 2 and 3 are all adiabatically tuned to zero in the same 
manner and the interqubit couphng between qubits 1 and 
2 is adiabatically switched on in the same way as that 
between qubits 2 and 3, then the left- and right-end Ma- 
joranas should occur with equal probabilities. Using this 
state = ^(7i^ + e*^7^)|F) as the initial state, one can 
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braid the left- and right-end Majoranas through the steps 
shown in Figs.[2fc)-|2fh) by adiabaticahy tuning the qubit 
parameters. This braiding of Majoranas corresponds to a 
unitary operator [3|,ll0| U = exp(7r7^7]^/4), which trans- 
forms -ff to U-ffU-^ = 7f and 7^ to U-f^U'^ = -jf. 
Therefore, the initial state of the system is trans- 
ferred to 1^)/ = |(7f - e'^-ff)\F) after braiding the 
left- and right-end Majoranas. 

Finally, we focus on probing Major ana fermions. The 
initial state = ^e'^j^)\F) given in Fig.Efb) is 
a ground state of the system with qubit V decoupled from 
the horizontal array of superconducting qubits. When ex- 
pressed in the basis states of qubits, this initial state can 
be written as = |^i23)i ^ where |^i23)i = 

AiilOiOsOs) + A,2|0i02l3) + A^alOihOs) + A^lOihla) + 
Ai5|li0203) + A,6|li02l3) + A,7|lil203) + Ai8|lil2l3), and 
= ^ii|Oi/) +^i2|li'). Also, the final state |^)/ = 
^{js ~ e*^7i^)|i^) is another degenerate ground state 
of the same system and can be expressed as |^)/ = 
1^123)/ ^ 1^1')/, where = ^/i|Oi/) + ^/2|li'). and 

1^123)/ has the same form as |^i23)i, but the A^^ are re- 
placed by Xfi^ with / = 1 to 8. The states and |^)/ 
can be distinguished using experimentally available state- 
tomography techniques for superconducting qubits (see, 
e.g., [ISl, ll9|), which involve reconstructing an unknown 
quantum state from a complete set of measurements of 
the system observables. 

Experimentally, it is more complicated to use state- 
tomography techniques to determine the quantum state 
of three qubits other than two qubits. Therefore, we 
can consider the state in Fig. [2fc) as the initial 
state. This state is a ground state of the system with 
qubits V and 1 decoupled from other qubits and can be 
decomposed as = |^23)i ^ ^ l^i)*^ where 

1^23) = Aa|0203) + A,2|02l3) + A,3|l203) + A,4|l2l3), 

=^ii|0i')+^i2|li'), and =7?ii|0i)+7?i2|li). 
From the state in Fig. [2fh), further proceeding with one 
step analogous to that from Fig. [2fb) to Fig. Efc), we 
achieve the final state with the originally unpaired Ma- 
joranas J2 and 7^ braided. This final state can also be 
decomposed as \^) f = 1^23)/ ^ l^i')/ ^ l^i)/^ where 
1*1'}/ = C/ilOi'} +e/2|li'), |*i>/ = r//i|Oi) +r?/2|li>, 
and 1^23)/ has the same form as |^23)i, but the Xu are 
replaced by A/^, with / = 1 to 4. Similarly, the states 
and 1^)/ can also be discriminated using state- 
tomography techniques. 

Discussion and conclusion. — When fabricating super- 
conducting circuits, parameter variations unavoidably 
occur, as in any solid-state system. For the charge-qubit 
array, Ej{^q) = can be achieved by having = ^^Oi 
irrespective of the parameter variations. Also, /i can be 
tuned, via the gate voltage Vg^ to the required value, even 
if Ech varies for different qubits. For varying Ejq among 
qubits, the inter qubit couplings also vary [see Eq. (|2j)]. 
One can replace the large Josephson junction by a dc 
SQUID and tune the SQUID, i.e., the effective Ejc, to 



obtain the desired value t for the interqubit coupling. As 
for the fiux-qubit array, e = can be achieved by hav- 
ing / = ^. Also, fi can be tuned to the given value 
by changing the fiux fs applied to the SQUID in each 
qubit. Moreover, even if the parameters of couplers vary, 
one can tune the flux fc in each coupler to achieve the 
required value t for the interqubit coupling [see Eq. ([7])]. 

In conclusion, we propose superconducting circuits to 
construct two superconducting-qubit arrays where Ma- 
jorana modes can occur. We suggest using four super- 
conducting qubits as the smallest system to demonstrate 
the braiding of Major ana modes, and show how to dis- 
tinguish the states before and after braiding Major ana 
fermions. These superconducting-qubit arrays can be ex- 
tended to wire networks, similar to the semiconducting 
wire networks in [lo|, to demonstrate the non-Abelian 
statistics of Majorana fermions. 
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